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Abstract 

Purpose- Landslides are major hazards to human activities, which often wreak havoc on economic resources, damaging 
properties and facilities in rural areas. The present study, considering that a perquisite of any development and planning is 
the recognition of the geographical features in an area, investigated the risk of landslide due to the expansion of agricultural 
land uses in rural areas. 

Design/Methodology/A pproach- This is an applied research that sought to examine the research background and select the 
most appropriate methods. Accordingly, it adopted a mixture of quantitative methods (fuzzy Delphi and fuzzy best-worst 
method), GIS and remote sensing techniques to achieve the research goal. 

Findings- According to the research findings, with increasing height, slope and vicinity to the fault lines, the risk of 
landslides rises in the study areas. These areas are mostly located in the highlands and the eastern and western regions, where 
rural areas are chiefly distributed. However, the majority of rural areas are distributed in the middle areas, which have better 
access to water resources and are in more favorable conditions due to topographic factors. Meanwhile, agricultural lands, 
due to the use of river water resources, have been distributed in the middle areas, which are classified as low risk areas in 
terms of landslides. In contrast, due to the limited flat lands in highlands, agricultural gardens have developed in highlands 
with a moderate slope, which subsequently pose the risk of landslide. Therefore, the regular monitoring of land use 
development to increase the safety factor in new housing construction and agricultural lands is one of the planning 
requirements for land use development in mountainous rural areas. 
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1. Introduction 
ccording to research conducted by the 
Center for Epidemiology, landslides 
account for about 17% of all deaths 
related to natural hazards worldwide 
(Kanungo et al., 2012; Pourghasemi et al., 
2012). Landslides are often recognized as 
one of the most devastating and widespread natural 
disasters in the world and a leading cause of death and 
economic losses (Achour et al., 2017; Chen et al., 
2018; Das et al., 2012; Pourghasemi et al., 2012). As 
Petley (2012) points out, landslide-associated losses 
and its destructive effects are more prominent in less 
developed countries, which is mainly rooted in the 
misconception of landslide hazards and the lack of 
appropriate resources (Razak & Mohamad, 2015). 
Landslides are more prevalent in mountainous areas 
and every year a host of landslides are reported in these 
areas (Lin et al., 2017; Shahabi et al., 2014). Given the 
future trends of these regions, as well as growing 
urbanization and unplanned = development, 
deforestation and growing regional rainfall due to 
climatic changes, especially in earthquake-prone areas, 
this natural phenomenon in expected to be intensified 
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in the future (Goetz et al., 2011; Kanungo et al., 2012). 
It would be of great concern, especially to mountain 
dwellers. 

As regards mountainous areas, it is obvious that the 
spatial deployment of rural areas in these areas is more 
affected by natural factors such as access to water, flat 
land and fertile soil, and therefore scant attention has 
been paid to factors such as natural hazards. As such, 
many rural areas are developed in areas that are 
vulnerable to natural disasters. In this regard, and 
shown by the research, most of the rural settlements in 
Tarom city are deployed in mountainous areas (Table 
1). In this context, during the past years, due to the 
prosperity of agricultural activities and the favorable 
environment of Tarom city for these activities, many 
agricultural activities and gardens have developed in 
valleys and the foothills. In some rural areas (Tahm, 
Chavarzagh, etc.), the sudden onset of landslides 
wreaks havoc on these agricultural uses. Accordingly, 
this study, recognizing the importance of land study 
prior to any planning for land use, investigates the 
development of agricultural land uses in rural areas of 
Tarom city with respect to the risk of landslides in 
these areas. 


Table 1. Distribution of residential areas (rural and urban) of Tarom city in natural types 


(Source: Statistic Center of Iran, 2011 & field observations, 2020) 


Area District County SUM Frequency Plain Frequency Mountainous Frequency 
Chavaisach Chavarzagh 30 %19.9 3 %10 27 %90 
Dastjerdeh 19 %12.6 6 %31.6 13 %68.4 
Rural Abbar 17 13 8 %AT A 9 %52.9 
Central Deram 41 %o27.2 =) %A2.2 36 %87.8 
Gilvan 42 %28.8 10 %23.8 32 %16.2 
Urban 2 %13 1 %o50 2 %50 
SUM 151 %100 33 %100 118 %100 

Apart from zoning areas for landslides, the lands to natural disasters, particularly in 


importance of addressing this issue is linked to the 
inadequate development of agricultural land uses, 
which are one of the main assets of rural 
households. 


2. Theoretical foundations of research 

Agricultural lands, which make up 40% of rural 
areas worldwide (Lesiv et al., 2019) are of the main 
land uses that have received growing attention due 
to the rising population and the need for food 
supply, development and exploitation of lands in 
rural areas. Meanwhile, despite the importance of 
land development and exploitation for agriculture in 
rural areas, it is necessary to recognize the use of 
land to mitigate the vulnerability of agricultural 
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mountainous areas that are prone to natural hazards. 
Landslides are one of the most prevalent natural 
disasters in mountainous areas. "Slide" is the motion 
of a mass under the impact of gravity, which is seen 
as a random process owing to the interaction of 
complex and unknown geographical, environmental 
and physical factors (Das et al., 2012). Landslides, 
unpredicted and destructive, are often considered as 
a natural hazard. Regarding landslides, risk 
assessment and zoning are the main measures in 
disaster risk management (Ambrosi et al., 2018) and 
one of the essential tools in any program. the main 
purpose of which is to mitigate the impact of natural 
disasters in the future (Skilodimou et al., 2019). 
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Over the past few decades, considering the 
importance of risk assessment of natural hazards 
such as landslides, a variety of methods have been 
developed for hazard mapping worldwide (Achour 
et al., 2017), among which GIS and remote sensing 
techniques have been widely utilized to assess areas 
more susceptible to landslide (Pirasteh & Li, 2017; 
Pirasteh et al., 2018; Shahabi et al., 2015). GIS 
(Kayastha et al., 2013; Wan et al., 2010), logistic 
regression models (Bui et al., 2016; Chen et al., 
2019; Das et al., 2010), bivariate and multivariate 
methods have been introduced in many studies as a 
suitable method for determining landslide 
susceptibility (Choi et al., 2012; Meinhardt et al., 
2015; Regmi et al., 2014; Zhang et al., 2016). 
Moreover, one way of evaluating a new approach is 
using the multivariate decision making methods, 
which is frequently used in research in combination 
with GIS method (Feizizadeh & Blaschke, 2014; 
Kayastha et al., 2013; Trinh et al., 2016). About 
Iran, this issue has been studied by Arab Ameri et 
al. (2018), Saffari & Hashemi (2017), Mansouri et 
al. (2016), Saffari (2014), and Moghimi et al. 
(2012). The Fuzzy Delphi decision making, AHP, 
ANP, Boolean logic and entropy methods have been 
employed to assess landslide risk. However, as the 
review of these studies suggests, they have 
primarily focused on zoning areas at the risk of 
landslide and scant attention has been paid to risk 
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Figure 1. Geographical location of the study area 
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assessment relative to land use development. With 
this in mind, the present study, by reviewing the 
research background and selecting the most 
appropriate methods, aimed at assessing the risk of 
landslide with respect to the development of 
agricultural land uses in rural areas. 


3. Research Methodology 

3.1 Geographical Scope of the Research 

The study area was Tarom County in Zanjan 
province, Iran. According to the latest political 
divisions of the country, this county comprises two 
central districts, Central and Chavarzagh. The 
Central districts comprises three villages of Abbar, 
Gilvan and Darram and Chavarzagh districts also 
includes two villages of Chavarzagh and 
Dastjerdeh. According to the last census (2016), this 
county had a population of 46641 people, of which 
219% lived in urban areas and 79% in rural districts 
of the county. In addition to the great geographical 
distribution of rural areas, the topographic type of 
rural and urban settlements in this city indicates that 
most settlements are located in mountainous and 
uneven areas. As regards employment and 
economic activities, considering the climatic and 
environmental conditions, horticulture, agriculture, 
services and industry sectors are the main source of 
employment, respectively. Figure 1 shows the 
geographical location of Tarom city. 
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3.2. Methodology 
The adoption of appropriate criteria for the 


intensity and susceptibility of landslides is a key 
step in its hazard analysis, on which the accuracy 
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of research results is dependent. Accordingly, in 
the first step of the research, the criteria affecting 
the risk of landslides are determined based on 
similar studies and climatic conditions of the study 
area, Table 2. 


Table 2. Criteria influencing landslides 


Factor Conditions of the study Area Background of the Research 
Height 
Slope Highlands in most of study areas (Blahut et al., 2010; Chen et al., 2017; Ghimire, 2011; 
Slope directions Kornejady et al., 2018; Zhang et al., 2015; Akgun, 2012; 
The Ghezel Ozan river running through Balteanu et al., 2010; Blahut et al., 2010; Chen et al., 
Distance from river the study area and natural river flowing 2017) 
in some regions 
ane i joenteta anual) area) Blahutet al., 2010; Akgun, 2012; Chen et al., 2017; Zhang 
earthquakes in 1990 ae) 
Geology Geological formation in most of studied (Balteanu et al., 2010; Chen et al., 2017; Ghimire, 2011; 
areas, low permeability to rainfall Zhang et al., 2015) 
Rain Rainfall and high relative humidity in the 
Relative humidity study area due to highlands and the (Balteanu et al., 2010; Chen et al., 2017; Zhang et al., 2015) 
Ghezel Ozan River 
lands Nios ands inch cha meeceaee (Balteanu et al., 2010; Blahut et al., 2010; Ghimire, 2011; 
barren and arid lands, Also, highlands in ogee etal) 
; Sou . (Blahut et al., 2010; Chen et al., 2017; Ghimire, 2011; Singh 
Vegetation the study areas are devoid of vegetation et al, 2014) 
Distance from road sisi aimee (Akgun, 2012; Chen et al., 2017) 


To screen the criteria, the experts’ opinion in fuzzy 
Delphi method was used. It was intended to draw on 
experts’ consensus regarding the selection of 
appropriate criteria as the research basis, and to analyze 
experts’ stances more precisely in a fuzzy space. The 
instrument was a researcher-made questionnaire based 
on fuzzy Delphi spectrum (Table 3). A total of 10 


university professors with a relevant educational 
background were randomly selected from the target 
population. The questionnaire first explained the main 
purposes of the research as well as the importance of 
accuracy in answering questions. Then, respondents 
were asked to state another factor related to the 
research goals besides the specified criteria. 


Table 3. Range of linguistic terms and numerical scale of fuzzy Delphi method 


Linguistic Terms Very Low Low Median High Very High 
ee (0, 0, 0.25) (0, 0.25, 0.5) (0.25, 0.5, 0.75) (0.5, 0.75, 1) (0.75, 1, 1) 


Appropriate screening criteria were selected and 
then based on the type of information obtained 
from each criterion, the zoning maps were drawn 
in ARC GIS, ENVI. These criteria were weighted 
using the best-worst method (BWM). This is one 
of the most effective methods for weighting the 
criteria, which was first proposed by Rezaei 
(2015). This method is superior to the AHP 
hierarchical method due to the compatibility ratio 
between the evaluation criteria. Given that this 
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method has a lower pairwise comparison, it 
provides more reliable results (Rezaei, 2015). In 
this method, first the most and least important 
criteria and sub-criteria in terms of the highest and 
lowest scores were determined using fuzzy Delphi 
method. Then, a _ pairwise comparison 
questionnaire was developed (Table 4). The 
specialized questionnaire was filled out by 15 
university professors acquainted with the research 
subject. 
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Table 4. Range of linguistic terms and numerical scale of the best-worst fuzzy method 


a Low importance Fair Very important Absolutely 

Dinenisietenied fee ance (WI) importance*(F1) (VI) important 
(EI) (AD 

Apw d,1,) (2.3, 1, 3.2) (3.2, 2, 5.2) (5.2, 3, 7.2) (7.2, 4, 9.2) 


After collecting the questionnaires, the data was 
fuzzily merged and its codes were weighed in the 
LINGO. After controlling the adjustment rate, the 
final weight of the criteria and sub-criteria was 
calculated. Then, by applying the weight of each 
criterion to the zoning map, the landslide risk map 
in the study area was drawn. However, to design a 
land use map for a 15-year period, Landsat 8 and 7 
images were obtained for the study area. After pre- 
processing the images and enhancing the spatial 


separation and mosaicization of the satellite 
images, the land use map was identified using the 
maximum probability classification method for 6 
land uses. After classification, for the verification 
and calculation of error matrix with terrestrial data, 
a comparison was drawn between samples of land 
use map and terrestrial data (GPS) (Figure 2). 
Finally, after overlapping landslide hazard map and 
land uses, residential areas and land use 
development were identified as landslide risk. 
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Figure 2. Land harvested points for validation of land use classification map 


4. Research Findings 
The target population was selected from among 
academic experts using fuzzy Delphi method. In 


n 
: 1 
aj = min{a;;} , bj = -) bij ‘ 
i=1 


Gg = max{ci;} , Crisp= 


this method, the questionnaire data was defined as 
triangular fuzzy numbers for each criterion and 
then the responses to the questionnaires were 
integrated according to Equation 1 
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Equation 1. Fuzzy Delphi method 


The median fuzzy value was set at 0.500 as the 
minimum fuzzy value to confirm the 
appropriateness of the criteria Accordingly, only 
the criterion of distance from main roads was 


removed as an inappropriate criterion (Table 5). 
Moreover, in the suggestions section, the experts 
approved the study criteria without offering any 
parameter. 
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Table 5. Fuzzy value of expert opinions and the appropriateness of each criterion 

Factor 5 4 eo ) 1 SUM Fuzzy value Review Consensus 
Height 1 9 - - [| - 10 (0.5, 0.77, 1) 0.757 Suitable %90 
Slope 1 7 2 - [| - 10 (0.25,0.71, 1) | 0.654 Suitable %T0 
Slope directions 2 8 - | - 10 (0.25, 0.54, 1) | 0.597 Suitable %T0 
Distance from fault 2 7 1 - | - 10 (0.25, 0.76, 1) | 0.671 Suitable %T0 
Rain 2 5 3 - | - 10 (0.25, 0.70, 1) | 0.651 Suitable %70 
Relative humidity - 3 7 - | - 10 (0.25, 0.56, 1) | 0.605 Suitable 50/: 
Distance from river 1 2 7 - - 10 (0.25, 0.58, 1) 0.610 Suitable 70/. 
Vegetation 2 2 6 - | - 10 (0.25, 0.62, 1) | 0.624 Suitable 60/. 
Geology - 4 6 - | - 10 (0.25, 0.59, 1) | 0.613 Suitable 60/. 
Land Use - 2 7 1 | - 10 (0,0.51, 1) 0.502 Suitable 70/. 
Distance from main roads - - - - 10 (0, 0, 0.75) 0.250 Unsuitable 40/. 

SUM 9 |) 43 | 51 | 5 - 7 . - 


The highlands are characterized with steep slope 
and instable foothills. Hence, in mountainous 
areas, mass displacements, especially landslides, 
are more probable. In terms of final weight, the 
physical criterion and land cover were determined 
as the most and the least important criterion, 
respectively. As regards physical criterion, the sub- 


criterion of height and slope directions were the 
most and the least important sub-criteria, 
respectively. As for hydrological criterion, rainfall 
and the distance from the river and as for land cover 
criterion, vegetation and land use were determined 
as the most and the least important sub-criteria for 
the landslide risk, respectively (Table 6). 


Table 6. Fuzzy Delphi final weight and determination of the most and the least important criteria by BWM 


method 
MethodBWM | Row | imal ao Sub-criterion pay al eae Factors 
Important 1 0.117 Height 
- 3 0.107 Slope ‘ 
Least important | _9 0,091 Sacco | One ya 
- 2 0.110 Distance from fault 
Important 4 0.107 Rain 
Least important 8 0.092 Relative humidity - 0.298 Hydrogeology 
- 6 0.100 Distance from river 
Important 5 0.102 Vegetation tas 
- 7 0.093 Geology important 0.278 Land cover 
Least important 10 0.082 Land Use 
- - 1 - - 1 SUM 


After identifying important and irrelevant criteria, 
the paired questionnaire was prepared by the BWM 
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method and filled out by 15 university professors, 
as shown in Figure 3. 
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Figure 3. Paired comparison of the best and worst alternatives using the BWM method 


From a gender perspective, the majority of experts professors. As for the field of study, most 
were male. They were chiefly in the age group of respondents were specialized in natural geography, 
40-60 years and held the position of assistant as shown in Table 7. 


Table 7. Demographic characteristics of the participants in the BWM 


Gender Frequency Percent Field of Education Frequency Percent 
Man 11 %73.3 Rural planning 3 %20 
Female 4 %26.7 Geomorphology 4 26.7% 
SUM 15 %100 Climatology 4 %26.7 
Age Group Frequency Percent Geology 4 %26.7 
More than 60 years 4 %26.7 SUM 15 %100 
40 to 60 years 7 %AG.7 Education Frequency Percent 
20 to 40 years 4 %26.7 Professor 2 % 13.3 
Less than 20 years 0 %O Associate 4 %26.7 
Not stated 0 %O Assistant Professor 6 %AO 
SUM 15 %100 PhD 3 %20 
SUM 15 %100 
Pairwise comparison judgments of criteria and sub- Then, according to Equation 2, in the LINGO 
criteria were merged in the form of fuzzy numbers, program, the fuzzy value of each criteria and sub- 
based on the most and the least important criteria is obtained. 


parameters, as shown in Table 8. 
Table 8. Fuzzy integration of criteria and sub-criteria based on Pairwise Comparison of the most and the least 
important criteria 


Factors Important Fuzzy integration Least important Fuzzy integration 
Physical d,1,) (0.67, 3, 4.5) 
Hydrogeology Physical (0.67, 2, 4.5) Land cover (0.67, 1, 2.5) 
Land Cover (0.67, 3, 4.5) d, 1,1) 
Sub-criteria Important Fuzzy integration Least important Fuzzy integration 
Height d,1,) (1.5, 3, 4.5) 
Slope : (0.67, 1, 2.5) bes (0.67, 2, 3.5) 
Slope directions Hage (15,3, 45) ee GLb 
Distance from fault (0.67, 1, 1.5) (0.67, 2, 3.5) 


23 


CN 
SRT? 


Journal of Research and Rural Planning No.2/ Serial No.33 
Factors Important Fuzzy integration Least important Fuzzy integration 
Rain (1,1, ) (0.67, 2, 4) 
Relative humidity Rain (0.67, 1, 2.5) Dp st (0.67, 1, 2.5) 
Distance from river (0.67, 2, 4) d, 1,1) 
Vegetation (1,1, ) (0.67, 1, 2.5) 
Geology Vegetation (0.67, 1, 2.5) Land Use (0.67, 1, 2.5) 
Land Use (0.67, 1, 2.5) (1, 1,1) 
mink min =k; 
(Ly, my, ul) ae [31 — 0.67 *uj3 = =k; lg, — 0.67 * uy = =—k; 
("mw") (15 j,Mpj,Upj)} S (k",k", k*) my, — 1 +m <=k;my, —1+m% > =-k; 
an a Sous Ug — 2.5 «15 S =k; ug, — 2.5 «15 = =—k; 
iy, my, uy — aw Mw, yw) S Cle", Ie", *) Iy, — 0.67 «uy < =k; ly, — 0.67 + uy > =—k; 
s.t. n my; 1* mj S =k;mg, —1*m3 2=—k; 
>. 8 (4) = ut, —2.5 «IY <=k;ug, —2.5 «2% > =—-k; 
= Iw — 0.67 «uy, < =k; 1 — 0.67 « ug, > =—k; 
iY < mj’ <u; mz — 1* M33 S =k;mj3 — 1* mg, = =—k; 
W >0 ul — 2.541% <=k;u —2.58+l% >=—-k; 
j=1,2,,n 0.167 « 1, + 0.668 * m, + 0.167 * u, + 0.167 « 


Equation 2. The best-worst method 


min =k; 

IB, — 0.67 * uj, < =k; lg, — 0.67 * uj, = =—k; 
mg, — 2* mj S =k; mpg, —2* mj = =—k; 

up, — 4.5*15 <=k;ug, —4.5 *]5 > =—k; 
If, — 0.67 «ujz < =k; lg, — 0.67 *ujz => =—k; 
mg, — 3 * mj < =k; mg, — 3 * m3 = =—k; 

upi — 4.5 * 3 S =k; ug, — 4.5 * 3 = =—k; 

| GS — 0.67 «ups < =k; US — 0.67 «ug, > =—k; 
my, — 1* mg3 < =k; my —1* mg; > =—k; 

uy, — 2.5 #05 <=k;uy —25*03 >=—k; 


0.167 « 1, + 0.668 * m, + 0.167 * u, + 0.167 * 


1, + 0.668 * mz + 0.167 * U> + 0.167 * 1; + 0.668 * 


m3 + 0.167 * u3; = 1; 
lL < my; m, < uj; 1, = 0; 
l, < Mz; Mz S Uy; ly = 0; 
l; < m3; m3 < uz; ly = 0; 


Equation 3. Main criteria in LINGO 
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I, + 0.668 * mz + 0.167 * U2 + 0.167 * 13 + 0.668 * 
Mz + 0.167 * uz = 1; 
lL, < my; m, < uy; Ll, = 0; 
l, = Mz; Mz = Uz; l, = 0; 
lz < m3; Mz S U3; 13 = 0; 


Equation 4. Land Cover sub-criteria in LINGO 


¥ — 0.67 «us <=k; 1%, — 0.67 « ull > 


mg,—1*m3< 
Up, — 2.5*13 < 
Ig, — 0.67 *uj3 = 
mpg, —2*m3S 
uy, — 418 < 
I) — 0.67 * up, < 
mj — 1*ms3 = 
<= 


w w 
ujs — 2.5 «Te 


ll 


min =k; 
=—k; 
k;mg, —1*m, ==—k; 
k 5 ug, —2.5*15 ==—k; 
k 51g, — 0.67 * uj >=-k; 
k;mg, —2*m3 2=—k; 
k;ug, —4*U3 ==—k; 
k ; Uj — 0.67 « ug, 2 =—k; 
k;mj3 —1* mg, = =—k; 
=k;ujz —2.5*]3 2=—k; 


0.167 « 1, + 0.668 « m, + 0.167 * u, + 0.167 « 
Il, + 0.668 * mz + 0.167 * uz + 0.167 *13 + 0.668 * 
M3 + 0.167 * uz = 1; 
lL, < my; m, < uy; Ll, = 0; 
l, < mz; Mz < Up; Iz = 0; 
lz < M3; M3 < Uz; Iz = 0; 


Equation 5. Hydrogeology sub-criteria in LINGO 
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min =k 
Ig, — 0.67 * ujs < =k; 1g, — 0.67 *uj3 => =—k 
mg, —1* my <=k;mZ, —1* mj => =—k 
ug, — 2.5 +05 <=k;ug, —2.5*U5 = =—k 
Ip, — 1.5 *ujg < =k 51g, — 1.5 * ujg = =—k; 
mg, — 3 *mjz < =k; mz, — 3 * mj = =—k; 
upi — 4.5 *]3 < =k; ug, — 4.5 «13 = =—k; 
IB, — 0.67 * ujg < =k 31, — 0.67 * uj, = =—k; 
ie Ay Sahm tee oa 
upi — 1.5 * U4 <=ks ug, — 1.514 => =—k; 
[jn — 0.67 * ups < =k; U7, — 0.67 « ugy = =—k; 
mjg — 2 * mB, < =k; mjz — 2 * mgs = =—k; 
uja — 3.5* U3 < =k; uj — 3.5 * 13 = =—k; 
[jg — 0.67 * ugg S =k; Ug — 0.67 * ugy = =—k; 
mj, — 2 * mes < =k;mjq —2*mz3 > =—k; 
uj, — 3.5 * 13 S=kj uj — 3.5 «13 => =—k; 


0.167 * 1, + 0.668 * m, + 0.167 * u, + 0.167 * 1, + 0.668 
* Mz + 0.167 * uz + 0.167 * 1, + 0.668 * mz + 0.167 * ug 
+0.167 * 1, + 0.668 * my + 0.167 * ug = 1; 

lL, < my; my, < uy; lL, = 0; 

l, < mz; Mz < Ug; Ip = 0; 
lz < mg; Mg < Ug; ly = 0; 
ly < m4; Mg < ug; ly = 0; 


Equation 6. Physical sub-criteria in LINGO 


The comparisons between criteria and sub-criteria our research method indicates its greater 
using the BWM are more consistent than other desirability. 
decision-making methods. The study of its value in 
Table 9. Fuzzy numerical value and compatibility of main criteria and sub-criteria 
Factors Fuzzy Weigh £ | CR | Row Sub-criterion Fuzzy Weight £ CR 
0250 [Slope (0147.0291,039) 
: ope 147, 0.291, 0. 
eset ee 3 Slope directions (0.082, 0.123, 0.209) DS | OSS 
; 4 Distance from fault (0.185, 0.291, 0.332) 
0.200 t t 5 Rain (0.344, 0.344, 0.757) 
Hydrogeology 0.225 S| Ss 6 Relative humidity (0.275, 0.275, 0.585) 0.068 | 0.068 
0.496 7 Distance from river (0.207, 0.207, 0.514) 
0.180 8 Vegetation (0.247, 0.356, 0.461) 
Land Cover 0.180 9 Geology (0.195, 0.356, 0.409) 0.026 | 0.059 
0.367 10 Land Use (0.174, 0.330, 0.330) 


Physical criteria had the highest weight for 


important and slope and land use had the least 


landslide risk zoning. Among the sub-criteria, significance for zoning the risk of landslides in the 
height and distance from the fault were the most regions (Table 10). 
Table 10. Standard weight and sub-criteria affecting the landslide risk 


Factors Weight Row Sub-criterion Weight Final Weight 
1 Height 0.314 0.164 
: 2 Slope 0.274 0.143 
Phys gmat 3 Slope directions 0.131 0.068 
4 Distance from fault 0.281 0.146 
5 Rain 0.415 0.111 
Hydrogeology 0.267 6 Relative humidity 0.327 0.087 
7 Distance from river 0.258 0.069 
8 Vegetation 0.356 0.075 
Land Cover 0.212 9 Geology 0.338 0.072 
10 Land Use 0.306 0.065 
SUM 1 - - - 1 
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For landslide riskoning, the final weight of the sub- 
criteria was applied to the zoning map of each and 
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areas exposed to landslide risk were identified 
(Figure 4). 


Legend 


® Rural Areas 
@ City Centers 
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In rural areas, a larger portion of land uses were 
dedicated to agriculture. Naturally, these land uses 
thrive in areas with a desirable access to water 
resources. Since the Ghezel Ozan river runs 
through the study area, most of the water areas are 
located in the middle areas. Residential areas, 
which are built in low and flat lands, and 


Figure 4. Landslide risk zoning 


re 


agricultural lands, due to the convenient access to 
river water resources, have been developed in the 
middle areas with a lower height and slope. 
However, the use of agricultural gardens (mostly 
olive) in the study period (2005-2009) was more 
prevalent in areas with medium and steep slopes 
(Table 11). 


Table 11. Land use area of areas classified by the maximum probability method 


Land Use 


Area (Hec) 


Changes 


ee) 


sum | amiss [ei [assis [soo |_| 


P= Kagem 88 | Kegncoecien re [= 
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Figure 5. Land use in the study area 


Given the overlap between the landslide risk and 
land use, the residential and agricultural land uses 
mainly developed in areas with a low landslide 


risk. The agricultural gardens primarily developed 
in areas with high risk of landslides, despite the fact 
that a landslide would damage these areas. 


Table 12. Land use development matrix in areas exposed to landslide risk over 2005-2019 


Landslide Area Residential Gardens of Agriculture Natural Barren 
ee = = sy ee ey a 


| Very Low __| Low 


| Median | * 


* 


a a 
a 9 ON 


5. Discussion and Conclusion 

Residents of rural areas in highlands often pursue a 
location-based biological pattern; however, due to 
the nature of their residence, they have to deal with 
a plethora of natural disasters in these areas. The 
rural areas under study sit in a mountainous area 
and landslides, as one of the natural hazards in rural 
areas, have always debilitated economic and 
infrastructural capacities of these areas. 
Considering that any development and planning 
requires knowledge and awareness of the 
geographical features of the region, which is a 
prerequisite for development, especially in rural 
areas, this study investigated the risk of landslides 
in light of the expansion of land uses. Given the 
purpose of the present study and as stated by 
Zumpano et al. (2018), this risk is largely ignored 
despite the fact that rural communities are 
vulnerable in terms of economic resilience and 
natural disaster losses. The present study adopted a 
combination of quantitative methods (fuzzy Delphi 
and best-worst fuzzy), GIS and remote sensing 
techniques to pursue the research goals by 
reviewing research background and selecting the 
most appropriate methods for zoning landslide risk. 
According to the results, physical criteria and sub- 


a ee ee 


criteria of height, distance from the fault and slope 
were the major landslide criteria for the risk of 
landslides (Arab Ameri et al., 2018; Saffari & 
Hashemi, 2017; Mansouri et al., 2016; Basharat et 
al, 2016). That is, with increased height, slope and 
proximity to fault lines, the risk of landslides 
surges (Skilodimou et al, 2018). These areas are 
located to the east and west of the study area. Since 
most of these areas are highlands, pastures and 
barren lands are the most common land. Moreover, 
the distribution of rural and residential areas in 
these areas is constrained and often temporary. In 
this context, considering that rural settlements 
usually develop around flat lands, with favorable 
soil and water sources, most rural areas have been 
distributed in the middle areas with a better access 
to river water resources (Ghezel Ozan) and most 
favorable topographic factors. Agricultural lands, 
due to the importance of access to river water, are 
distributed in the middle areas, which run a lower 
risk of landslide. In contrast, agricultural gardens, 
due to constraints related to flat and even lands in 
rural areas, have usually expanded in high and 
medium-slope areas. However, given that the 
analysis and identification of areas for any 
development and planning is one of the key steps 
to hamper financial losses on land use in the event 
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of hazards, it is important to allocate a greater 
attention to this issue. Therefore, some practical 
ways are suggested to reduce losses and financial 
losses associated with landslides for agricultural 
gardens in rural areas: 
eConsidering the economic value _ of 
agricultural lands and proper development of 
these land uses in areas with the risk of 
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e Regular monitoring of land use development 
to enhance safety in new housing construction 
and agricultural lands. 

e Providing the necessary infrastructure to raise 
awareness of rural residents about the dangers 
of landslides. 

e Developing rural infrastructure services in 
low and flat areas to diminish landslide 


landslide, 
e Evaluating and identifying arable lands in 
areas with a lower risk of landslides. 


susceptibility. 
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